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ABSTRACT 

We present results from a pilot HST ACS deep imaging study in broad-band V of five low- 
redshift QSO host galaxies classified in the literature as ellipticals. The aim of our study is to 
determine whether these early-type hosts formed at high redshift and have since evolved passively, 
or whether they have undergone relatively recent mergers that may be related to the triggering of 
the nuclear activity. We perform two-dimensional modeling of the light distributions to analyze 
the host galaxies' morphology. We find that, while each host galaxy is reasonably well fitted by a 
de Vaucouleurs profile, the majority of them (4/5) reveal significant fine structure such as shells 
and tidal tails. These structures contribute between ~5% and 10% to the total y-band luminosity 
of each host galaxy within a region of r ~ 3 r e g and are indicative of merger events that occurred 
between a few hundred Myr and a Gyr ago. These timescales are comparable to starburst ages in 
the QSO hosts previously inferred from Keck spectroscopy. Our results thus support a consistent 
scenario in which most of the QSO host galaxies suffered mergers with accompanying starbursts 
that likely also triggered the QSO activity in some way, but we are also left with considerable 
uncertainty on physical mechanisms that might have delayed this triggering for several hundred 
Myr after the merger. 

Subject headings: galaxies: active - galaxies: interactions — galaxies: evolution — quasars: general 
1. INTRODUCTION 
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Quasi-stellar objects (QSOs) are the most lu- 
minous active galactic nuclei (AGNs), believed to 
be powered by accreting supermassive black holes 
(BH). The growing observational evidence that 
massive BHs exist not only in the centers of active 
galaxies, but also in inactive early-type galaxies 
indicates that more than the mere presence of a 
massive BH is needed to trigger the activity. And 
while the "bright quasar phase" is now often cited 
as an essential phase in th e evolution of galaxies 



(e.g.. lSpringel et alJl2005al >. we still do not have a 



clear picture of what causes the triggering, or even 
what the necessary conditions are for the trigger- 
ing to occur. 
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For more than two decades now, mergers have 
been suspe cted to be re s ponsible for trigger - 
ing QSOs (|Stocktonl Hjjii Zanders et all Il988h . 
There is an abundance of circumstantial evi- 
dence in the literature connecting QSO activity 
at lo w redshifts with mergers and interactions 
Heckman et al.l[l98l Irlutchings et alll 19881 



Canalizo fc Stockton! I200I lUrrutia et al.1 I2OO7F 
Not only do a number of low-redshift QSOs show 
signs of interactions and mergers, but there is also 
a tendency for QSO host galaxies t o exist in only 



moderately rich environments (e.g., [ Dressier et al- 
ii 98. 4 lEllingson et al.l[l99~D: IWold et alJl200ll see. 



however, Martini et al.1 2007 for a counterexam 
pie) consistent with the results of quasar cluster- 
ing studies, which suggests that q uasars exist in 
haloe s of mass ~10 12 - 10" M ()Hopkins et al 
20071 Ida Angela et all 120071 ). In such environ- 



ments, velocity dispersions are moderate, and in- 
teractions have their greatest effect. 

At least in the case of QSOs found in ultra- 
luminous infrared galaxies (ULIRGs), there is a 
close connection between mergers and QSO activ- 
ity. ULIRGs are essentially always fou nd in merg- 
ing systems ( Sanders fc Mirabel 19961 ). and far too 
many of them harbor QSOs to be explained by the 
joint probability of finding a given L* or brighter 
galaxy to be both a ULIRG and a QSO h ost galaxy 
by chance (jCanalizo fc Stockton! l200lh . So, at 
least for this subclass of QSO, some of the gas 
that the merger feeds into the central starburst 
also finds its way much deeper to the central black 
hole. 

Nevertheless, recent high-resolution imaging 
studies show that a large fraction of QSOs reside 
in hosts with r elaxed morphologie s like elliptical 



galax i es (e.g.. iDisnev et al.l 11995 



19971 : iMcLure et al.l 119991 : iFlovd et all 



Bahcall et al.. 
20041 ). 



These observations have elicited claims that tidal 
interactions in QSO host galaxies are the excep- 
tion rather than the rule, and that even in those 
cases in which there are clear signs of interac- 
tion, the interactions are no t necessarily related 
to the nuclear activity (e.g., iDunlop et al.l I2003L 
hereafter D03). From two-dimensional fitting of a 
well matched sample of 33 AGNs (radio-loud and 
radio-quiet QSOs as well as radio galaxies), using 
images obtained with the Wide Field & Planetary 
Camera 2 (WFPC2) on board HST, D03 con- 
cluded that "for nuclear luminosities My < —23.5, 



the hosts of both radio-loud and radio-quiet AGN 
are virtually all massive elliptical galaxies with 
basic properties that are indistinguishable from 
those of quiescent, evolved, low-redshift ellipticals 
of comparable mass" . 

Such conclusions conflict with recent findings 
of significant amounts of dust and cold molec- 
ular g well as young stellar populations, 
in many rad i o-loud and radio-quiet QSO s (e.g . 
Martel etal.1 Il999l: ICanalizo fc Stockton! 12 001 



Evans et al.l200ll:IScoville et al.l2003l:lTadhunter et al 
20051 : iBarthell 12006^ . Although the hosts of 
most luminous AGNs are bulge-dominated, they 
seem to be significantly bluer than inactive el- 
liptical galaxies and indeed show evidence for 
starbursts in the relatively rec e nt past (1-2 
Gyr) (e.g., iKauffmann et alJl2003l : I Sanchez et al.l 



2004; Canali zo et al.ll2006l : iSchweitzer et al.ll200' 
Jahnke et al.ll2007t ). 

In addition, both numerical simulations and 
observations show that merger remnants that 
have undergone violent relaxation have surface- 
brightness profiles that follow r 1 / 4 laws (e.g. 



Toomre fc Toomrdl97ilBarnesll98ilHibbard fc van Gorkuml 
19961 ). even while sti ll showing clear signs o f 
the tidal interactions (|Rothberg fc Joseph! 120041 ). 
Thus, if QSO host galaxies are indeed relatively 
recent merger remnants, they would most likely 
have elliptical profiles, but they would also show 
fine structure indicative of past interactions. 

Detecting fine structure in the already elusive 
QSO host galaxies is a challenging task. Studies of 
nearby merger remnants show that their fine struc- 
ture is often faint compared to the galaxies them- 
selves. For example, the total luminosity of the 
shells in shell galaxies accounts only fo r 5-15% of 
the t o tal lum i nosity of the galaxy (e.; 



19861 : IPrieurl 119881 IWilkinson eta! 



.Fort et al 
20001 ) 



and 

their detection often requires a resolution of a t 
least ^0.5 kpc (e.g., Schweizer fc Seitzer 1992). 
Therefore, detecting fine structure in QSO hosts 
requires high angular resolution and higher signal- 
to-noise (S/N) observations than those obtained in 
previous studies. 

Thus, we are conducting a study with deep (5 
orbit) HST ACS and WFPC2 observations of a 
sample of QSO host galaxies that are classified 
as ellipti cals. We presented r esults for the first 
object in lCanalizo et al.1 (|2007l . hereafter Paper I). 
In the present paper, we describe results for the 
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four remaining objects from our pilot ACS study 
of 5 objects. We will present results for the rest of 
the sample (14 additional objects) in subsequent 
papers. Throughout this paper, we adopt fl\ = 
0.7, Q m = 0.3 and H = 71 kms" 1 Mpc" 1 . 

2. SAMPLE SELECTION 

Since the goal of our project is to investigate 
the possibility that QSO host galaxies are formed 
through mergers even when they do not show dra- 
matic signs of tidal interactions, we should ideally 
draw our sample from a sample of hosts that have 
been classified as undisturbed elliptical galaxies. 
The best such sample is that of D03. 

From the detailed morphological study of radio- 
loud and radio-quiet QSO host galaxies by D03, 
we selected QSOs for which the host galaxies were 
classified as elliptical galaxiefl Moreover, we se- 
lected those objects for which existing deep Keck 
spectra did not reveal extended emission lines that 
could contaminate our observations in the broad 
T^-band filter (F606W) . We chose a sample of five 
QSOs, three radio-quiet and two radio-loud. The 
sample properties are summarized in Table Q] 

3. OBSERVATIONS AND DATA RE- 
DUCTION 

All observations were obtained using the ACS 
Wide Field Channel (WFC) onboard the HST 
with the broad F-band F606W filter (AA = 
2342A; 1 pixel corresponds to 0.05"). Five sets 
of dithered images were taken (with four subsets 
of 485-588 s integration time each), yielding a to- 
tal integration time of 10740-11464 s per target 
(see TableQ}. 

In order to achieve background-limited images 
enabling us to detect any existing faint structure 
in the QSO host galaxies, we obtained individual 
frames with long integration times. Consequently, 
the central pixels of the QSO nuclei are saturated 
in these images. We discuss the effects of the sat- 
uration in detail in the next section. 

We re-calibrated the data manually, correct- 
ing the pipeline flat-fielded single exposures for 



the bias-level offset between the adjacent quad- 
rants that is still present in the fi nal product 



1 Note that the D03 radio- loud and radio-quiet QSOs were 
selected to be statistically indistinguishable in their optical 
luminosity and redshift. 



Pavlovskv ct al 



Koekemoer et al 



2005). We used 



2002) to combine 



of CALACS 
MultiDrizzle 
the individual images to the final distortion and 
cosmic-ray corrected scientific image (default val- 
ues plus bits=8578 and a deltashift file with the 
offsets between the images as determined from 
stars within the field-of-view (FOV)). Note that no 
sky subtraction was performed during MultiDriz- 
zle. 

4. IMAGE ANALYSIS 

Different methods suc h as unsharp masking , 
creating a structure map (jPogge fc Martinill2002Q . 



and using the t wo-dimensional ga laxy fitting pro- 
gram GALFIT (jPeng et alj|2002h were applied to 
look for any fine structure that might be present. 
We here describe only the procedure we adopted 
for the use of GALFIT as the other two methods 
did not yield additional information (see Paper I 
for details on these methods). 

For convolution with the point-spread function 
(PSF) of the HST ACS optics, we created both an 
artificial PSF star from TinyTinU (Version 6.3) at 
the same position as each QSO and a PSF star 
from an observed star on an ACS/ WFC F606W 
image observed close to the positions of the QSOs. 
For the Tiny Tim PSF, we adopted a power law 
F v = v a with a = —0.3 which is the average value 
we obtain for the QSOs[f| The computed radius 
for this artificial PSF is 10". 

For the observed PSF star, we chose a star 
with a signal-to-noise ratio (S/N) of 20,000 from 
the HST archive. The star was observed on 
20 dithered images with a total exposure time 
of -8100 s (GO-9433, data sets j6mfl9* and 
j6mf21*) and at a position of the WFC chip less 
than 5" away from the position of each QSO; note 
that the full ACS FOV is 3.4' x 3.4'. These images 
were processed in the same manner as described 
above for the QSO observations. Finally, we elim- 
inated a few faint objects surrounding the PSF 
star and modified the PSF image to minimize any 
noise introduced during PSF subtraction and con- 



Available at 

http: / /www. stsci.edu/software/tinytim/tinytim. html 
3 (-B — R) ranges between 0.1 and 0.7 as measured from nu- 
clear Keck spectra. 
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Table 1 
Details of Observation 



QSO 

(1) 


IAU designation a (J2000) 
(hh mm ss) 
(2) (3) 


S (J2000) 

(° ' ") 
(4) 


z 

(5) 


Type 
(6) 


Exp. Time 

(sec) 
(7) 


scale 
(kpc/") 
(8) 


PHL909 


0054+144 


00 57 09.9 


+ 14 46 10 


0.172 


RQQ 


10920 


2.885 


PKS 0736+01 


0736+017 


07 39 18.0 


+01 37 05 


0.191 


RLQ 


11364 


3.137 


PG 0923+201 


0923+201 


09 25 54.7 


+ 19 54 05 


0.190 


RQQ 


11464 


3.124 


MC2 1635+119 


1635+119 


16 37 46.5 


+ 11 49 49 


0.146 


RQQ 


11432 


2.520 


OX 169 


2141 + 175 


21 43 35.5 


+ 17 43 49 


0.211 


RLQ 


10740 


3.392 


Note.— Col. (1): 


QSO. Col. (2): 


IAU designation. 


Col. (3,4): 


Optical positions 


taken from NED 


. Col. (5): 



Heliocentric redshift as listed in NED. Col. (6): Classification as either radio-quiet QSO (RQQ) or radio-loud 
QSO (RLQ). Col. (7): Total integration time in sec. Col. (8): Physical scale in kpc/ . 



volution procedures. To adaptively smooth the 
PSF, we compared data values to the standard 
deviation s of the surrounding sky: (1) for data 
values > 7s, we retained the unmodified PSF; (2) 
for data values between 3s and 7s, the image was 
smoothed with a Gaussian kernel with a = 0.5 
pixel; (3) for data values < 3s, a Gaussian ker- 
nel of a = 2.0 pixel was used, and (4) for data 
values < Is after this last operation, the value 
was replaced with 0. The radius of the observed 
PSF star is ~4". The full-width-at-half-maximum 
(FWHM) of the PSF is 0.1" in both the artificial 
star and the observed PSF star, in close agreement 
with the FWHM of stars on the ACS images close 
to the QSO. 

The observed PSF star has (B — R) ~ 1.4 mag, 
as inferred from USNC0 data. To test chromatic 
effects on the form of the PSF, we created a series 
of artificial PSF stars using TinyTim. We used 
power laws with a ranging from -4.0 to 0.9, corre- 
sponding to (B — R) values between 2 mag and 
mag (to account for the observed range of (B — R) 
of the QSOs and the star, plus possible uncer- 
tainties). We also created PSFs using a black- 
body spectrum with temperatures between 3000 
and 40000 K. We fitted these different PSFs and 
a host galaxy model to our QSOs. The derived 
host galaxy parameters vary by 2%. We thus con- 
sider the chromatic effects as negligible compared 
to the other uncertainties such as the PSF mis- 
match due to telescope or plate scale breathing, 
focus changes, PSF saturation (see below) etc. 

We subtracted the two different PSFs (star and 
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TinyTim) from saturated and unsaturated stars 
within the ACS FOV, for comparison. The ob- 
served PSF star yielded significantly better results 
than the TinyTim PSF. Thus, we used the ob- 
served PSF star for all the subsequent fitting and, 
in the following, we simply refer to it as PSF or 
PSF star. This exercise also allowed us to estimate 
the quality of our PSF fit, taking into account the 
saturation of the QSO nucleus. The central region 
within a radius of ~1.7" is strongly affected by the 
PSF subtraction. Any structure within this region 
may be an artifact. 

We created artificial images to test the influ- 
ence of the PSF saturation on the derived QSO 
host galaxy parameters. Since luminous quasars 
generally do not vary much over time scales of a 
few years, we can compare our quasar luminosi- 
ties with D03 to estimate the degree of saturation 
in our images and the potential systematic errors 
that might arise due to saturation. This direct 
comparison shows that our AGNs may be satu- 
rated by as much as 60%. (Note that we define 
"60% saturation" in such a way that 60% of the 
flux is lost due to saturation.) In addition, we cre- 
ated model images where the AGN has the same 
luminosity as in D03. We then saturated the AGN 
by 60%. After masking out the central saturated 
pixels in the same manner as we analyzed our data 
(see below) , we find that our systematic errors can 
be as large as 7% for the magnitude, 30% for the 
effective radius of the different host galaxy compo- 
nents and up to 30% for the Sersic index. The PSF 
magnitude can be underestimated by up to 90% 
due to saturation. This explains the differences in 
AGN luminosity that we observe compared with 
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D03 (see Ta ble Hj). 

GALFIT (jPeng et alJl2002h is a 2-dimensional 



fitting program based on x 2 minimization; it al- 
lows the user to simultaneously fit one or more 
objects in an image with di fferent mode l light dis- 
tributions J j3uch_as_Sersic JSersig ll968[) , de Vau- 
couleurs ( de Vaucouleurs! 19481 ). or exponential. 
The following steps were performed. We created a 
mask to exclude the saturated pixels in the center, 
the diffraction spikes, surrounding objects, and 
any visible fine structure in order to fit only the 
smooth underlying host galaxy light distribution. 
For the purpose of the x 2 minimization, pixels 
were weighted by 1/c 2 with a based on Poisson 
statistics as estimated by GALFIT. (Note that the 
weighting is however irrelevant for the central sat- 
urated pixel that were masked.) A PSF and a 
Sersic function were fitted. We fitted the back- 
ground sky and bright close neighboring galax- 
ies simultaneously, while faint neighboring objects 
were simply masked. This least-square fit was then 
subtracted from the original image to obtain the 
residual image in which any fine structure super- 
posed on the smooth host galaxy light distribution 
is more readily visible. 

To estimate the uncertainties caused by the sky 
determination, we used two different approaches. 
First, different sizes of the fitting region were used 
(up to the largest possible size) with the sky vary- 
ing freely. Second, the sky was also determined 
independently and held fixed to this value during 
alternative fits. 

Note that a Sersic power law is defined as 



S(r) = Seffexp 



where £ e ff is the pixel surface brightness at the ef- 
fective radius r e s, and n is the Sersic index (n = 4 
for a de Vaucouleurs profile, n = 1 for an expo- 
nential profile). 

To fit the host galaxies, we used three differ- 
ent models: (1) a single de Vaucouleurs profile 
(Sersic with n = 4), (2) a single Sersic profile (n 
free) and (3) a combination of a de Vaucouleurs 
plus exponential (Sersic with n = 1) profile. If 
the sky value is kept fixed during the fitting pro- 
cess, the fit for (1) has a total of nine free pa- 
rameters: three parameters for the fitting of the 
PSF (a, 5, magnitude), and six free parameters 




for the fitting of the de Vaucouleurs profile (a, 5, 
magnitude, r c g, b/a, P.A.). Model (2) has one 
additional free parameter for the index n. Model 
(3) has six more free parameters than model (1). 
These parameters correspond to the exponential 
profile (a, 5, magnitude, r g, b/a, P.A.). In all 
three cases, if the sky is also fitted, there are three 
additional parameters (background and gradient 
in a and S). In some cases, we used a Sersic pro- 
file to fit bright galaxies in close proximity to the 
QSO hosts (PHL909: six galaxies; PG 0923+201: 
seven galaxies, MC2 1635+119: one galaxy). In 
those cases, the free parameters for these fits add 
to the number of total free parameters. 

5. HOST GALAXY PROPERTIES AND 
FINE STRUCTURE 

In this section, we summarize the results for 
all five QSOs in the pilot sample. Although 
MC2 1635+119 was discussed in detail in Paper 
I, Figure [3] shows its luminosity profile which was 
not included in Paper I. To facilitate comparisons, 
we also include the results for MC2 1635+119 in 
the tables. 

5.1. Host Galaxy Morphology 

Figure Q] shows the ACS/WFC images for 
four of the five QSOs in our sample: PHL909, 
PKS 0736+01, PG 0923+201, and 0X169. In- 
cluding MC2 1635+119, four of the five host galax- 
ies reveal extended fine structure on different 
scales and with different morphologies, ranging 
from spectacular shells and arcs (MC2 1635+119; 
Paper I), to tidal tails (PHL909) to jet- like 
(OX 169) and spiral-like (PKS 0736+01) structure. 
Only one object (PG 0923+201) resides in a host 
galaxy without any obvious fine structure. We 
here describe the overall host-galaxy morphology 
for each object in turn. Further details on each 
object, including a comparison with results in the 
literature and a description of neighboring objects, 
are given in Appendix [A"l 

The host galaxy of PHL909 shows a variety 
of fine structures at different radii. A central 
ring- or disk-like structure (r ~1.8" - 2.5") can 
be seen, possibly intersected by a dust lane to 
the SE (Fig. [I]). To the NW, two shells are ap- 
parent. They occur at distances of 2.5" (P.A. ~ 
328) and 3" (P.A. ~ 298). Diffuse outer material 
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seems to form another ring or disk-like structure 
(Fig. H inset; r ~ 3.4" - 6"). The most promi- 
nent fine structure consists of two tidal tails, one 
extending from the disk-like structure in the SE, 
curving towards the N, and the other extending 
from the disk in the NW, bending towards the 
south (Figs. Q]and|4|). The western tidal tail ex- 
tends to a galaxy ("a"; Fig. g} ~16" to the NW 
(P. A. ~ 284°). The tidal feature most likely corre- 
sponds to the faint e xtended emission reported by 
iMcLure et al .1(1 19991) which was a lso seen in the K- 
band image ( Dunlop et al.lfl993T ). The companion 
("a") itself has two tidal tails, one reaching out in 
the direction of the QSO and connecting with its 
tidal tail, another tidal tail extending to the S up 
to 5.5" (Fig. [T]). Also the eastern tidal tail seems 
to be connected to a very faint small companion 
("d") at 14" distance from the QSO (P. A. -98°), 
possibly a dwarf galaxy (Fig. |4j inset). 

PKS 0736+01 has been observe d in th e near- 
infrared K-band by iDunlop et all (|l993h . The 
host is about 10" in diameter and is described as 
highly disturbed with a large area of low surface- 
brightness nebulosity and several compact com- 
panions embedded. In our deep HST ACS im- 
age, this nebulosity is clearly visible and shows re- 
solved structure that looks like a faint face-on spi- 
ral. The spiral-like structure extends out to a ra- 
dius of 16" (~ 50kpc), i.e. very large for a normal 
spiral galaxy. Also, the pitch angle changes with 
radius, arguing against a simple spiral-arm inter- 
pretation. Instead, what is observed here, may be 
spatial wrapping of material from a merger event 
(see discussion in Sec. [6]). 

PG 0923+201 resides in a galaxy group of at 
least 6 galaxies, some of which form closely in- 
teracting pairs, including the two large elliptical 
galaxies to the south-east of the QSO (see Ap- 
pendix [Aj. However, the host galaxy itself does 
not reveal any obvious fine structure in the deep 
HST ACS image. 

The host galaxy of MC2 1635+119 has been de- 
scribed in detail in Paper I. It reveals spectacular 
fine structure, consisting of (at least) five inter- 
leaved inner shells (out to r ~ 13 kpc) and several 
extended arcs out to a radius of ~ 65 kpc. 

The host galaxy of OX 169 exhibits an extended 
linear str ucture which has been reported by several 
authors (IStocktonl Il978l: Irlutchings et all 1 1984J: 
Gehren et alJl984tlSmTth et al.ll986UrIeckman et al 



19861) . OX 16 9 was studied bvlStockton fc Farnham 



(|1991I ) and iHutchings fc Nefil dl99J ) who de- 



Stockton & Farnham 



scribe d this jet-like structure in detail 
(|1991l ) showed that it consists mainly of old stars 
at the same redshift. This is supported by the 
pro minence of the struct ure also in the K-band im- 
age ( Dunlop et al.lll993h . While the feature resem- 
bles a galactic disk seen edge-on, with a dust lane 
in the SE, it is not centered on the QSO; almost 
2/3 of it lies on the SE sid e of the QSO. Because 
of its strong asymmetry, IStockton fc Farnham 
( 199ll ) favor an interpretation as a tidal tail seen 
edge-on rather than a normal galactic disk. In 
our deep HST ACS image, the linear structure 
is very prominent and extends over a projected 
distance of ~68kpc. The feature runs from NW 
to SE (P.A. ~ 303°), out to -7.6" in the NW 
and ~12.5" in the SE. It is slightly curved and 
appears to continue across the nucleus. A dust 
lane is clearly visible between 7.6" and 11.5" in 
the SE side of the tail. In addition to the promi- 
nent linear structure, the deep HST ACS image 
reveals faint extended structure for the first time: 
a shell-like feature in the NE (P.A. ~ 83°) at a 
distance of ~5", plumes to the N of the SE end of 
the tail (P.A. ~ 110°) at a distance of ~11.5" from 
the center, and a fan extending from the SE tail 
towards the SW. The latter extends out to a com- 
panion ("a"; Fig. HJ) to the SW (P.A. ~ 235°) at a 
distance of 6.1" and may be a tidal stream. The 
companion has an overall spiral-like morphology 
which seems to be tidally stretched. These find- 
ings provide strong support for the interpretation 
of the linear structure as a tidal tail seen edge on 
rather than it being a normal galactic disk. It is 
only natural to assume that all features have a 
common origin during a merger event. 
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Fig. 1.— ACS/WFC images of four QSO host galaxies, all smoothed with a a = 0.5 pixel Gaussian. In 
each image, north is up and east is to the left. A harder stretch is shown as inset, Gaussian smoothed with 
(7=1 pixel (PHL909, 0X169) and a = 2 pixels (PKS 0736+01, PG 0923+201), respectively. Including 
MC2 1635+119 (Paper I), fine structure can be seen in at least four of the five QSO host galaxies (PHL909: 
tidal tails and shells; PKS 0736+201: spiral-like structure; 0X169: disk-like structure, most likely tidal tail 
seen edge-on; MC2 1635+119: inner shells and outer arcs, see Paper I). 
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5.2. Host Galaxy Properties 

As described in Section EJ we modeled the 
QSO host galaxies with GALFIT and fitted sin- 
gle de Vaucouleurs profiles, single Sersic compo- 
nents, and a combination of de Vaucouleurs and 
exponential-disk profiles. Varying the range of the 
sky parameter shows that the uncertainty in pho- 
tometry can be up to 20%, the uncertainty in r e g 
up to 10%, and the uncertainty in Sersic index 
up to 20%. These uncertainties add to the un- 
certainties from the PSF saturation (see Sec. HJ. 
The total uncertainties are ~20% for photometry, 
~30% in r c ff and ~35% in the Sersic index. The 
results are summarized in Table [5] 

For all objects, good fits were achieved with 
all three fitting functions mentioned above. The 
fits improved slightly in going from a single de 
Vaucouleurs to a single Sersic to a combination 
of de Vaucouleurs plus exponential-disk profiles 
(see Table [2]). While this is expected simply be- 
cause of the increasing number of free parame- 
ters, the residual images do look smoother for the 
two-component fits. In such a decomposition, the 
de Vaucouleurs profile remains the dominant com- 
ponent, but the exponential-disk component con- 
tributes between 6% and 28% of the total luminos- 
ity in V of the hosts (Table [2]). The residual image 
did improve significantly when adding a second 
component to the host in MC2 1635+119 (Paper 
I)- 

Single-Sersic fits yielded a Sersic index around 4 
for each of the host galaxies (ranging between 2.9 
and 6.5). Our results for the surface-brightness 
profiles of the underlying host gal axies are thus 



i n ove rall agreement with those of iDunlop et al 
(|2003l ). In Figure [2l we show the residual images 



derived by subtracting a GALFIT model consist- 
ing of a single Sersic profile from the observed im- 
ages. 



FIT and its subcomponents. The resulting ID 
profiles for the different fits (single de Vaucouleurs, 
single Sersic, de Vaucouleurs plus exponential) 
are only marginally different. Therefore, we only 
present the best fitting model using a Sersic com- 
ponent in Figure [3] 

The deviation of the model from the data 
within the central 1.7" is most likely caused by the 
PSF subtraction. This is especially pronounced 
for PG 0923+201, as this object has the brightest 
PSF (see Table [2}, and the low S/N ratio in the 
outer parts of the PSF star results in additional 
noise in the data. In general, data points at a 
higher surface brightness than the model corre- 
spond to the fine structure that was masked out 
during GALFIT fitting but that was included for 
the creation of the ID profiles. This is clearly the 
case for OX 169, for which the tidal tail causes the 
deviation between data and model beyond 5". 

The surface-brightness plot of the observed host 
profile of PG 0923+201 shows a bump peaking 
around 2" that is not fitted by the model. Since 
this feature lies just outside of the central region 
that is affected by the PSF, we found it difficult 
to determine whether it is real or not. However, 
it could reflect the presence of an inner disk. If 
so, the outer radius of this disk would be approx- 
imately 8kpc. The outer "uplift" of points over 
the best-fit profile (around r ~ 8") is most likely 
due to the companion galaxy to the SSE and the 
relatively strong residuals that remain after sub- 
tracting the GALFIT model for this galaxy (see 
also Fig. [2]). The observed surface brightness for 
PHL 909 also deviates slightly from the fit with a 
bump peaking around 2.3". This bump coincides 
with the ring-like structure seen in the image, with 
an outer radius of ^7kpc. 



5.3. Surface-Brightness Profiles 

We derived surface-brightness profiles for each 
object using the IRAI0 task ellipse. The pro- 
files were obtained for both the observed surface 
brightness and the best-fitting model using GAL- 



5 IRAF (Image Reduction and Analysis Facilities) is dis- 
tributed by the National Optical Astronomy Observato- 
ries, which are operated by AURA, Inc., under cooperative 
agreement with the National Science Foundation. 
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Table 2 

Results of Modeling the QSO Host Galaxies Using GALFIT 





ivioclcl 


Function 


Z-i (Ct,0 ) 

(", ") 


"1F606W 


1 I 


feff 






Scrsic 


U In 
V j II 


r.J\. 


_2 

x„ 


(fitted area) 






(mag) 


(mag) 


(") 


(kpc) 




index 




(°) 




(1) 

V- 1 -/ 


V W 


\''f 


(4) 


(5) 


(6) 


(71 


(8) 




(9) 


(10) 


(11) 


fl2l 


PHL909 


deV+Exp 


PSF 


(0,0) 


17.23 
















4.83 


(950x950) 




deV 


(-0.1,-0.1) 


17.08 




2.35 


6.77 


4 


(fixed) 


0.57 


130.4 








Exp 


(-0.96,1.45) 


19.88 




3.14 


9.05 


1 


(fixed) 


0.57 


111.0 






S 


PSF 


(0,0) 


18.18 
















5.14 






S 


(0.1,-0.06) 


16.79 




2.53 


7.30 




5.6 


0.63 


130.2 






dcV 


PSF 


(0,0) 


17.39 
















5.17 






deV 


(-0.1,-0.05) 


16.96 




2.54 


7.32 


4 


(fixed) 


0.61 


130.4 






dcV (D03) 


PSF 






15.5 




















deV 






16.6 


2.76 


7.93 


4 


(fixed) 


0.61 


131 






S (D03) 


S 














4 








PKS 0736+01 


deV+Exp 


PSF 


(0,0) 


17.18 
















2.20 


(1051x1051) 




deV 


(0.01,-0.05) 


17.70 




1.63 


5.10 


4 


(fixed) 


0.94 


2.3 








Exp 


(-0.94,-0.15) 


18.75 




7.09 


22.23 


1 


(fixed) 


0.80 


86.2 






S 


PSF 


(0,0) 


16.07 
















2.25 






S 


(-0.04,0) 


17.08 




4.27 


13.41 




6.5 


0.98 


32.7 






doV 


PSF 


(0,0) 


16.86 
















2.27 






deV 


(-0.04,0) 


17.35 




3.30 


10.34 


4 


(fixed) 


0.98 


39.2 






doV (D03) 


PSF 






16.2 




















deV 






16.9 


3.27 


10.21 


4 


(fixed) 


0.97 


26 






S (D03) 


S 














5.3 








PG 0923+201 


deV+Exp 


PSF 


(0,0) 


16.39 
















2.59 


(1051x1051) 




deV 


(-0.04,-0.02) 


17.40 




1.44 


4.49 


4 


(fixed) 


0.96 


174.2 








Exp 


(1.02,0.75) 


20.45 




2.26 


7.07 


1 


(fixed) 


0.51 


146.1 






S 


PSF 


(0,0) 


16.15 
















2.76 






S 


(0.07,-0.03) 


17.48 




1.81 


5.66 




2.9 


0.97 


42.1 






doV 


PSF 


(0,0) 


16.36 
















2.78 






deV 


(-0.06,0.01) 


17.31 




1.68 


5.2 


4 


(fixed) 


0.97 


42.9 






dcV (D03) 


PSF 






15.7 


















deV 






17.2 


2.02 


6.30 


4 


(fixed) 


0.98 


38 






S (D03) 


S 














3.3 








MC2 1635+119 


deV+Exp 


PSF 


(0,0) 


17.58 
















3.49 


(481x481) 




deV 


(-0.03,0.04) 


17.50 




1.12 


2.82 


4 


(fixed) 


0.71 


57.9 








Exp 


(0.53,0.36) 


18.81 




5.48 


13.8 


1 


(fixed) 


0.82 


27.7 






S 


PSF 


(0,0) 


17.02 
















3.91 






S 


(-0.03,0.04) 


17.18 




2.47 


6.23 




4.95 


0.71 


55.9 






doV 


PSF 


(0,0) 


17.17 
















3.93 






dcV 


(-0.03,0.04) 


17.27 




2.05 


5.16 


4 


(fixed) 


0.71 


56.2 






dcV (D03) 


PSF 






18.1 




















deV 






16.8 


2.28 


5.73 


4 


(fixed) 


0.69 


56 






S (D03) 


S 














5.6 









Note. — Col. (1): QSO and, in parentheses, the size of the fitted area in pixels. Col. (2) GALFIT model (dcV — dc Vaucoulcurs, Exp 
— Exponential, S — Sersic). Col. (3): Individual components used. Col. (4): Offsets with respect to the PSF. Col. (5): Integrated 
apparent magnitude in the F606W filter. Note that the magnitudes mF606W given for the PSF underestimate the AGN component as 
the PSF was saturated in all cases; see text for details. Col. (6): Integrated apparent magnitude converted from F675W to Cousins R 
band as given by D03. Col. (7,8): Effective radius in " and kpc, respectively. Note that although we concluded in Section \4\ that the 
central region within a radius of ~1.7" is still strongly affected by the PSF subtraction, all r G ff given in this table are more than 10 
times the full- width- at- half- maximum (FWHM) of the ACS PSF. Col. (9): Scrsic index. Col. (10): Axis ratio. Col. (11): Position angle 
(east of north). Col. (12): \ 2 P cr degree of freedom, v. \t corresponds to the model in column 2. 

Results from D03 are listed for comparison. The r e ff given here for the D03 results were recalculated according to the worldmodel 
adopted in this paper. The P. A. given here for the D03 results were derived by adding the P. A. given in their Table 3 to the orientation 
of the spacecraft (as the P.A. given in their Table 3 is apparently not corrected for it). 
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Table 2 - continued 



QSO 


Model 


Function 


A(a,8) 
(, ") 


?™F606W 


R 




r e ff 


Sersic 


b/a 


P.A. 


xl 


(fitted area) 






(mag) 


(mag) 


(") 


(kpc) 


index 




(°) 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


OX 169 


doV+Exp 


PSF 


(0,0) 


17.80 














3.45 


(900x900) 




deV 


(0,0) 


17.68 




1.82 


6.17 


4 (fixed) 


0.47 


120.4 








Exp 


(0,0.15) 


19.34 




1.23 


4.17 


1 (fixed) 


0.66 


18.5 






S 


PSF 


(0,0) 


17.54 














3.66 






S 


(0,0) 


17.42 




1.30 


4.41 


3.9 


0.64 


122.9 






doV 


PSF 


(0,0) 


17.56 














3.66 






deV 


(0,0) 


17.41 




1.29 


4.36 


4 (fixed) 


0.64 


123.0 






dcV (D03) 


PSF 






15.9 


















deV 






17.3 


1.88 


6.35 


4 (fixed) 


0.47 


121 






S (D03) 


S 












3.6 









Note. — Col. (1): QSO and, in parentheses, the size of the fitted area in pixels. Col. (2) GALFIT model (deV — dc 
Vaucoulcurs, Exp — Exponential, S — Sersic). Col. (3): Individual components used. Col. (4): Offsets with respect to the 
PSF. Col. (5): Integrated apparent magnitude in the F606W filter. Note that the magnitudes mF606W given for the PSF 
underestimate the AGN component as the PSF was saturated in all cases; see text for details. Col. (6): Integrated apparent 
magnitude converted from F675W to Cousins R band as given by D03. Col. (7,8): Effective radius in ' and kpc, respectively. 
Note that although we concluded in Section^] that the central region within a radius of ~1.7 ' is still strongly affected by the 
PSF subtraction, all r c ff given in this table are more than 10 times the full-width-at-half-maximum (FWHM) of the ACS PSF. 
Col. (9): Sersic index. Col. (10): Axis ratio. Col. (11): Position angle (east of north). Col. (12): \ P er degree of freedom, v. 
X„ corresponds to the model in column 2. 

Results from D03 arc listed for comparison. The r c ff given here for the D03 results were recalculated according to the worldmodcl 
adopted in this paper. The P.A. given here for the D03 results were derived by adding the P.A. given in their Table 3 to the 
orientation of the spacecraft (as the P.A. given in their Table 3 is apparently not corrected for it). 



Table 3 

Contribution of Fine Structure to Host Galaxy Light 



QSO 


fs (3r cff ) 


S fs (3r e e) 


Ef s (max) 




% 


mag arcsec -2 


mag arcsec -2 


(1) 


(2) 


(3) 


(4) 


PHL 909 


7.8 


29.2 


24.8 


PKS 0736+01 


8.2 


29.6 


26.5 


MC2 1635+119 


6.3 


28.3 


25.2 


0X169 


13.3 


29.5 


23.7 



Note.— Col (1): QSO. Col. (2): Percentage of flux con- 
tained in the observed fine structure within 3r e ff- Col. 
(3): Surface- brightness magnitude of the fine structure in the 
F606W filter within 3 r c ff . Col. (4): Average surface-brightness 
magnitude within an area of 2" x 2" around the maximum of the 
fine structure. 
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Fig. 2. — Residual images derived by subtracting a GALFIT model consisting of a single Sersic profile. In 
each image, north is up, east is to the left. MC2 1635+119 was presented in Paper I. 
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Fig. 3. — Surface-brightness plots of observed and best-fitting model profiles for all five QSO host galaxies. 
The observed profile is shown as diamonds with error bars. The resulting profile of fitting a single Sersic 
profile is shown as a solid line. The residuals (fit — data) are shown in the lower panels. 
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5.4. Fine Structure 

To estimate the luminosity within the vari- 
ous types of fine structure compared to the to- 
tal luminosity of the host galaxy, we created a 
mask that included all the light within an annu- 
lus with an inner radius 1.7" and an outer radius 
3r e ff (as taken from a single de Vaucouleurs fit, 
see Table [5]). The image was multiplied by this 
mask (good=l, bad=0), and the total counts in 
the product were summed. This was done both 
for the GALFIT residual image (/a ne structure) ob- 
tained by subtracting the GALFIT model of a 
de Vaucouleurs + exponential profile, and for the 
GALFIT model itself (/ ga i axy ). We then computed 
the ratio /finest ructurc / /galaxy The fractional lumi- 
nosity of the extended fine structure in such an 
annulus ranges between 2% and 5%. To account 
for regions that were over-subtracted by GALFIT, 
we also used a slightly different approach: Apply- 
ing the same general method as above, we only 
summed the pixels that are > — la of the sky back- 
ground. Then, the fractional luminosity of the ex- 
tended fine structure ranges between 6% and 13%. 
The results are summarized in Table |31 

However, note that these values give only the 
average fraction within 3 r c ff . Locally, the surface- 
brightness magnitude of the fine structure rela- 
tive to the underlying model can be much higher. 
The surface-brightness magnitude within an area 
of 2"x2" centered on the maximum of the fine 
structure lies between ^24 and ~26 magarcsec -2 
(as derived from the GALFIT residual image) for 
the four objects with fine structure (Table [3]). 

6. DISCUSSION 

We find significant fine structure in the deep 
HST ACS images of at least four of the five 
QSO host galaxies (PHL909, PKS 0736+01, 
MC2 1635+119, and 0X169) in our pilot sam- 
ple. While the prominent linear structure in 
OX 169 h as been previously reported by several 
authors dStocktonl Il978t iHutchings et al.l 11984 



Gehren et al.lll984lsin"ith et al.ll986tlHeckman et al 
1986t iMcLure et alj I1999T ) , the structure in the 
other three objects is revealed here for the 
first time. Only one of the objects studied, 
PG 0923+201, does not show any obvious fine 
structure. 

Our findings are supported by the WFPC2 im- 



ages of D03 and IMcLure et all (|l999h since some 
of the structure that we find is also visible in their 
images: Part of the north- western tidal tail i n 
PHL909 is visible (Fig. All. lMcLure et al.lll999T ): 
some of the spiral-like s tructure in PKS 073 6+01 
can be seen (Fig. A6 in McLure et al. 1999); and 
the shell in MC 1635+119 labeled "d" in Fig. 3 
of Paper I is visible to the north-east of the 
QSO nucleus at a distance of ~4" (Fig. A18 in 



McLure et al. 1999f ). However, the depth reached 



in our ACS images is needed to positively identify 
the "fuzz" seen in these earlier images as real. Our 
results clearly show the need of high S/N imaging 
for an accurate interpretation of QSO host galaxy 
morphologies. 

The host galaxy of PG 0923+201 does not show 
obvious signs of fine structure. However, the 
surface-brightness profile reveals a bump in the 
inner 2" that could reflect the presence of an in- 
ner disk. Interestingly, PG 0923+201 resides in 
a galaxy group of at least 6 galaxies. It is the 
only QSO in our sample that has been observed 
with the Spitzer Space Telescope, both with the 
Infrared Spectrograph (IRS) and the Multiband 
Imaging Photometer for Spitzer (MIPS; 24/im). 
Neither FIR emission nor Polycyclic- Aromatic 
Hydrocarbon (PAH) features have been detected 
dShi et al.l l2007h. interestingly, this lack of signs of 
star-formation activity in the infrared goes hand- 
in-hand with the lack of obvious signs for recent 
merger activity. 

Much of the fine structure observed in our 
deep HST ACS images of the other four QSO 
host galaxies is indicative of merger events in 
the relatively recent past. The possible type and 
timescales for the merger in MC2 1635+119 were 
discussed in detail in Paper I. Briefly, we found 
that the fine structure seen in the host galaxy of 
this object was most likely produced by a strong 
tidal interaction within the last ~1 Gyr. On the 
one hand, the regularity of the inner shell struc- 
ture is suggestive of a near radial collision of a 
dwarf galaxy with a giant elliptical galaxy in a mi- 
nor merger. On the other hand, the arcs observed 
at much larger distances than the inner shells, and 
other tidal debris off-axis from the direction of the 
encounter implied by the inner shells, are more 
indicative of a major merger. Using simple N- 
body simulations, we estimated that the time scale 
needed to form the observed structure is less than 
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~1.7 Gyr (see Paper I for details). In compari- 
son, the spectrum of the host galaxy is dominated 
by a population of intermediate-age (~1.4 Gyr) 
stars, indicative of a strong starburst that may 
have occurred during the merger event. Thus, the 
observed QSO activity may have been triggered 
in the recent past either by a minor merger or by 
debris from an older (~Gyr) major merger that is 
currently "raining" back into the central regions 
of the merger remnant. 

The prominent extended structure observed in 
0X169 is most likely a tidal tail seen edge on. 
Prominent tidal tails are also observed in the 
host of PHL909. The tails in both QSO hosts 
are strongly suggestive of relatively recent major 
merger events, where at least some of the parent 
galaxies were disk galaxies. 

The structure surrounding PKS 0736+201 re- 
sembles that of the Seyfert 1.5 galaxy NGC5548, 
an Sa galaxy with ripples, tidal arms and a faint 
long tidal tail that has bee n described in detail by 
Schweizer fc Seitzerl ([1988). The observed struc- 
ture in NGC 5548 has been interpreted in terms of 
the spatial wrapping of material from an obliquely 
infalling companion. When seen in projection, 
such material can app ear in the form of re l ativel y 
sharp-edged ripples (jSchweizer fc Seitzerl Il988| ). 
Compared with NGC 5548, the faint spiral-like 
structure seen in PKS 0736+01 is more extended 
(~ 50kpc vs. lOkpc) and more regular. 

Structure similar to that observed in PKS 
0736+201 can be reproduced in numerical simula- 
tions of minor mergers. For example, the merger 
of a dwarf elliptical galaxy with a large spiral 
and a total mass ratio of 1:8 results in very sim- 
ilar structure ~1 Gyr after the first passage of 
the two galaxies, but also again at lat er stages in 
the merger event ([Younger et alJl2007t T. J. Cox, 
private communication). 

By comparing the structure observed in these 
hosts to results from numerous published iV-body 
simulations of mergers (e.g. 



Hibbard fc Mihol 



19951), we know that the ages of these tidal fea- 
tures are likely of the order of a few orbital times. 
Given the masses of giant elliptical galaxies and 
the distances at which we observe these features, 
this time span ranges from a few hundred Myr to 
about 1 Gyr. In the cases in which we can clearly 
see tidal tails, we can infer that the dynamical 
timescales cannot be much less than ~ 250 Myr, 



since that is roughly the amount of time that it 
would take for stars with typical orbital velocities 
(~ 200 kins -1 ) to travel in a straight line to the 
projected maximum distances at which we observe 
the features (e.g., - 50 kpc for PHL909 or ~ 42 
kpc for 0X169). On timescales longer than ~1 
Gyr, tidal tails are expected to rapidly disperse 
and lose contrast. 

Further evidence that these are relatively 
evolved mergers is the fact that enough time has 
elapsed since the initial tidal encounter for each 
host to acquire an essentially elliptical morphol- 
ogy. While it is possible that some of the ob- 
jects, such as MC 1635+119 and PG 0923+201, 
may have started off as elliptical galaxies, at least 
0X169 and PHL909 must have resulted from the 
merger of disk galaxies. The fact that they are 
at the present time reasonably well fitted by de 
Vaucouleurs profiles implies merger timescales on 
the order of a Gyr. 

To address the question of whether the AGN 
activity may indeed have been triggered by the 
merger events, we need to obtain estimates of the 
different timescales involved. How do the ages 
of the stellar populations compare to the time 
elapsed since the merger events, and what does 
that imply for the QSO duty cycle? What is the 
connection between merger, star formation, and 
BH accretion? 

As mentioned above, the timescales for the 
merger events observed in our sample is of the 
order of a Gyr or less. G. Canalizo & A. Stock- 
ton (2008, in preparation) find similar timescales 
for strong starburst events in the host galaxies of 
these QSOs. By modeling deep Keck spectra of 
these objects, they find traces of major starburst 
episodes in intermediate-age starburst populations 
that involve a substantial fraction of the stellar 
mass. These starburst ages range from ^0.7 Gyr 
(0X169) to 2.2 Gyr (PKS 0736+017). Hence, it 
seems possible that the starbursts were triggered 
by the merger events. 

The most difficult timescale to estimate is that 
of the duration of the AGN activity. Current es- 
timates of QSO lifetimes range between 10 6 and 
10 8 years and are deri ved primarily from demo- 
graphics (see review bv lMartinill2004). For exam- 
ple, fr om theoretical calculations Yu fc Tremainei 
(2002) estimate a mean lifetime for a luminous 



QSO (iboi > lO^ergs" 1 ) of (3-13) x 10 7 



yr 
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which is comparable to the Salpeter time scale (~ 
4.5 x 10 7 yr for e == 0.1 and L/L e dd = l: ISalpeter 
1964 : Martinil 2004 ). Apart from demographics, 



there have been several attempts to deduce the 
AGN duty cycle for individual galaxies. For ex- 
ample, the length of radio jets can yield a lower es- 
timate of QSO lifetimes, de pending on the expan- 



sion speed of the jet (e.g., Willis fc Strom 1978; 



IScheuerlll993lBlundell et al.lll999l ). However, this 



method cannot be applied to our sample as the two 
radio-loud QSOs in the sample each feature only 
a compact flat-spectrum radio source. Also, the 
size of the NLR has been proposed as a str aightfor- 
ward geometric measure of AGN lifetimes (jMartinil 
120041 ). However, the NLR sizes for the five QSOs 
of our sample have not yet been determined. 

Thus, as we do not have a direct way to esti- 
mate the time elapsed since the triggering of the 
nuclear activity for the objects in our sample, it is 
difficult to work out the physical relationship be- 
tween the mergers we detect and the fueling of the 
QSOs. An overly simplistic scenario, in which each 
QSO was triggered by the merger at the time of 
the starburst episode would imply QSO lifetimes 
on the order of a Gyr, in stark contrast to the the- 
oretical estimates cited above. Instead, our obser- 
vations may provide evidence in support of signifi- 
cant time delays between tidal interactions and the 
fueling of the central black holes, as predicted b 



hydrodynamic simulations (see, e.g.. 



Springel et al. 2005bl ; Hopkins et al 



Barnes 199 



2007). 



However, while QSOs often reside in bulge- 
dominated galaxies like the ones presented in this 
study, they are also frequently found in young 
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mergers (e.g., Urrutia et al.l|2007|: iGuvon. Sanders k, Stockton! 



2006 



ICanalizo fc Stockto n 2001; Hutching s et al 



19941 ). If, as the accumulating evidence suggests, 
mergers are indeed essential for the triggering of 
QSOs, then our results imply at least one of two 
scenarios: Either QSO activity is episodic (e.g., 



Norman fc Scovillel[l988l;lHopkins et al.ll2006h and 
occurs over longer timescales than previously spec- 
ulated, or there is a large range of values for the 
time delays between the merger and the onset of 
the nuclear activity. 
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A. NOTES ON INDIVIDUAL OBJECTS 



Here, we give a more detailed description of each of the four QSOs. (Note that MC2 1635+119 was 
discussed in detail in Paper I.) Specifically, we review and discuss published data, compare them with our 
findings, and include a description of neighboring objects. The redshifts of companion galaxies and other 
galaxies within the ACS FOV are taken from NED unless stated otherwise. The full ACS/WFC images are 
shown in Figure SI with labels for the galaxies to which we refer in the following discussion. Note that we 
cannot show detailed images for each such galaxy; instead, we refer the interested reader to the HST archive. 



A.l. PHL909 



Based on older images, the host g alaxy of the radio-quiet QSO PHL 909 (z=0 . 172, 1" ~ 2.89 kpc 



was described as an elliptical galaxy ( Bahcall et al. 1996L 1997 : Taylor et al. 19961 : McLure et al. 199 



>0 
)9: 



Hamilton et al.l I2002T ) . The age of the stellar popula tion from off-nucl ear optical spectra is not well con- 



strain ed, but the spectra are suggestive of an old age (jNolan et al.ll200lh . Note, however, that lHughes et al 



(2000) report a significant increase of flux bluewards of a weak 4000-A break, but cannot determine whether 
it is due to scattered light from the QSO or indicative of a younger stellar population. 

BH mass estimates (making use of a, FWHMha and Ismn, i?-band l uminosity of the host spheroid) cove r 
a range of 0.1 7 -2.5 x 10 9 M ? (jMcLure fc Dunlo^feoOlllWu et alJliooHDunlop et1iil|2003l : IWu fc Liulbooi . 
Dunlop et all (|2003l ) predict an Eddington ratio (L nU c/iedd) of 5.2%. IChun et all (j2006T ) classify PHL909 
as Lyman limit system. They mention a number of faint small objects surrounding the Q SO that could 
be co mpanions. In a K-band image, the galaxy extends ~10" towards a western com panion (Dunlop ct al. 
1993). The optical spectrum of PHL909 shows doub le-peaked Balmer e mission lines ( Strateva et al.l 20031 ). 



1998). 



In the radio, PHL 909 is an unresolved point source (jKukula et al. 

The arc- like features at ~2" from the QSO mentioned by Bahcall et al. (1996) correspond most likely 
to the centra l ring- or disk-like structure (r ~1.8" - 2.5") seen in our deep HST ACS image. However, 



iBahcall et al. ( 19961 ) were unable to determine whether the features were artifacts due to the much lower 
S/N of their images. 

The neighborhood of PHL 909 is populated with several objects, some of which look like dwarf galaxies. 
Seven brighter, nearby galaxies make PHL 909 appear to reside in a galaxy group. However, redshifts are 
known for only 3 objects, of which only one object is at the same redshift as the QSO, the large spiral galaxy 
("h") to the S (P. A. ~ 183°) 28" away (z = 0.169). The two other objects are foreground galaxies, one an 
elliptical galaxy ("b"; z = 0.102) to the NE (PA. ~ 28°) at a distance of 18", and one a spiral ("e"; z = 
0.124) to the SE (PA. ~ 140°) 20" away. Two lenticular objects with unknown redshift lie at 31" E ("c"; 
PA. ~ 96°) and 12" S ("P; PA. ~ 173°) projected distance of the QSO. A spiral galaxy to the SW ("g"; 
PA. ~ 200°) is 20" away. Another bright object is the galaxy that is apparently interacting with the QSO 
(see Section IBTTj) . In the full FOV of the ACS frame (Fig. [4j), redshifts are known for two other galaxies, both 
foreground spirals ("i", z = 0.104, PA. ~ 29.3°, distance ~ 1.12' and "j", z = 0.082, PA. ~ 76°, distance 
~ 2.34'). 



A.2. PKS 0736+01 

PKS 0736+01 (z=0.191, 1" ~ 3. 14 k pc) was classified as a blazar bv lAngel fc Stockman! (jl980h . as a 
compact flat - spectr um radio quasar bv IWall fc Peacock! (|1985l ). and a s an optically violen t ly variable by 
Brown et al.l (19891). Th e host galax y was described as an elliptical bv lWright et al. ( 1998 ): McLure et al 



(|l999|) : lFalomo fc U lrich (2000), and Hamilto n et al 



(12002 ). Off- nuclear optical spectra have been obtained 



200lh . indicating an age of 12 Gyr. The estimated BH 



and fitted bv lHughes et all (|2000l ) and lNolan et all 
mass (from relations wi th a, FWHMhh and £51 nn 

0.14 and 1.3 x 10 9 M (|McLure fc Dunlopll200lUWu et alJliool iDunlop et al.ll2003l ). IDunlop et al 



.R-band luminosity of the host s pheroid) lies between 



give 3.5% as an estimate of the Eddington ratio. The optical HST WFPC2 image of iMcLure et al 



2003) 



1999) 



shows three companions (two to the NW and one to the SE), but fails to reveal the low surface-brightness 
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nebulosity, described as highly disturb ed, seen in the K-band ( Dunlop et al. 19931 ). In the near- infrared 



H-band, the host is resolved but round ( Kotilainen et al. 1998) 



Several faint objects surrounding PKS 0736+01 underlie the isophotes of the "spiral" . Two of these are 
background objects as determined from Keck spectra (G. Canalizo & A. Stockton 2008, in preparation): a 
lenticular-shaped galaxy at 5.9" to the SE ("a"; Fig.[4] P.A. ~ 147°; z = 0.8500) and an irregular galaxy at 
a distance of 4.4" to the NW ("b" ; P.A. ~ 307°; z = 0.3785). The latter shows a tidal tail that is extending 
toward another object to the NW ("c"; P.A. ~ 297°; distance ~ 9") of unknown redshift. 

Whether the observed faint spiral-like structure in PKS 0736+01 indeed indicates the presence of a faint 
spiral host galaxy, whether it is formed from debris of a merger event as speculated in Section O or whether 
it is actually a combination of both (a spiral disturbed by accreted material), our observations clearly show 
that the host galaxy of PKS 0736+01 cannot be considered to be an undisturbed elliptical galaxy. 

A.3. PG 0923+201 

The radio-quiet QSO PG 0923+201 (z=0.19, 1" ~ 3.12kpc) resides in an elliptical host ga laxy 
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(|Bahcall et al.lll997t McLure et al.lll999Ll2000| lHamilton et al.ll2002t iGuvon. Sanders fc Stocktonll2006l). It 



i s not detected as a radio source (jKukula et al.lll998l ). From a noisy optical off- nuclear spectrum. iNolan et al 



(|200ll ) estimate a stellar-population age of about 12 Gyr. BH mass estim ates (using a, FWHMha and L 5W q 
E-band luminosity of the h ost spheroid) yield M RH ^ 0.11-2.6 x 10 9 M Q (jMcLure fe Dunlodl200 ltlWu et al 



20021 iDunlop et alJliool . iDunlop et all(|2003l ) estimate an Eddington ratio L nuc /L cdd of 8.8%. 

PG 0923+201 is a member of a galaxy group : In t he AC S FOV, at least 5 gala xies are at a comparable 



redshift of z ~ 0.19, (pair "a" . iHeckman et al 



Ellingson et all (|199lh : pair "d" 



Keck spectra; 
"d"; Fig. UK. 



(|1984T ): "c" 

G. Canalizo & A. Stockton 2008, in preparation), four of which form two interacting pairs ("a" 
The closest pair ("a") at ~ 12" to the east consists of two large elliptical gala xies [one galaxy at 10" (P.A. ~ 
153°) and one galaxy at 15" (P.A. ~ 161°) to the SE] . iHutchings et al.l(|l989h speculated that the QSO may 
be interacting with these. The northern galaxy of this pair shows a shell-like structure. The other pair ("d") 
lies at a projected distance of ~ 1.6' to the north-east (P.A. ~ 21°). Both galaxies have clear tidal tails on 
one side and are connected by a tidal bridge. The fifth galaxy ("c") at a comparable redshift lies 18" to the 
east (P.A. ~ 83°), a small galaxy with overall elliptical morphology but faint extended tail-like structures. 
The tidal tails make it resemble a more distant and faint version of the Antennae galaxies (NGC4038/4039), 
the best-studied and nearest (~ 19.2 Mpc) example of a major interaction between two massive gas-rich 
spiral galaxies. However, this object is possibly at a later stage of the merger as the two nuclei seem to have 
coalesced. The SE arm ends in a knot, a compact object that could be an interacting dwarf galaxy. 

Note that within the ACS FOV, another galaxy pair ("e") at a distance of ~ 2' to the east (P.A. ~ 92°) 
is at an advanced stage of merging, wit h curved tidal tails and the nuclei appearing disrupted in several 
fragments. While lEllingson et al. ( 199ll ) list this pair at a redshift of z = 0.1899, Keck spectra show that 
it is rather at a redshift of z = 0.2318 (G. Canalizo & A. Stockton 2008, in preparation). The redshift is 
known for one other galaxy within the ACS FOV, a background spiral ("f" , z = 0.2328, P.A. ~ 87°, distance 
~ 2.6')- There is another very close galaxy to the SW (P.A. ~ 224°) of the QSO, ("b"; distance ~ 11") with 
peculiar, tidally disrupted morphology, but unknown redshift. 

There are at least two arcs close to the QSO (labeled as "1" and "2" in Fig. [4]) that could be gravitationally 
lensed galaxies. However, they may also be tidal debris, especially in the case of "1" , which lies close to the 
interacting pair marked "a" . 



A.4. OX 169 

The h ost galaxy of the radio-loud QSO 0X169 (z=0.211, 1" ~ : 3.39 k pc) has been descri bed as being an 
elliptical (jHutchings fe Nefll992UTavlor et al.lll996l : lMcLure et al.lll999l ). INolan et alj (|200lh infer an old age 
for the stellar population of OX 169 from their optical off-nuclear spectrum, although the spectrum has a low 
S/N ratio. BH mass estimations using different relat ions (a, FWHMha and Lsmn, i?-band luminosity of the 
host spheroid) cover a range of 0.22-1.7 x 10 9 M Q (jMcLure fc DunlopllioOlt IWu et alll2002l : IDunlop et al ' 



17 



20031 ). The Eddington ratio is estimated to be 6.1% dDunlop et al.ll2003l). From its radio emission, OX 169 



is classified as a compact (< 4") flat-spectrum source ( Feigelson et al.l[l984l ). 

0X169 seems to have several faint satellite galaxies, but their redshifts are unknown. The closest one at 
3.7" to the NE ("c"; P.A. ~ 43°) merges with the outer isophotes of the QSO host galaxy and resembles a 
dwarf elliptical galaxy. Two others lie towards the SW at distances of 4.6" ("b"; P.A. ~ 245°) and 6.1" ("a"; 
P.A. ~ 235°). The former has an irregular shape and resembles a tidally disrupted galaxy, while the latter 
object is the companion already described in Section f5.ll 
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Fig. 4. — Full ACS/WFC images of four QSO host galaxies. The two bars mark the position of the QSO. 
North is up, east is to the left. The FOV of the ACS/WFC frame is 3.4' x 3.4'. As inset, a smaller area 
around each QSO is shown, with the image size identical to the images in Figs. [1] and O Letters mark the 
galaxies we refer to in the text. The letters refer to the closest object north or south of them. In the case of 
PC 0923+201, the two numbers refer to the arcs east of them. 
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